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ABSTRACT 

We report the detection of the eclipse of the very-hot Jupiter WASP- 12b 
via z'-band time-series photometry obtained with the 3.5-meter ARC telescope 
at Apache Point Observatory. We measure a decrease in flux of 0.082±0.015% 
during the passage of the planet behind the star. That planetary flux is equally 
well reproduced by atmospheric models with and without extra absorbers, and 
blackbody models with / > 0.585±0.080. It is therefore necessary to measure the 
planet at other wavelengths to further constrain its atmospheric properties. The 
eclipse appears centered at phase = 0.5100+o:oo6L 

consistent with an orbital 

eccentricity of |ecosci;| = 0.0161;q;qq9 (see note at end of §4). If the orbit of the 
planet is indeed eccentric, the large radius of WASP- 12b can be explained by 
tidal heating. 
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Introduction 



The transiting Very Hot Jupiter WASP-12b, discovered by lHebb et al.l (|2009[ ). has many 
notable characteristics. With a mass of 1.41 ± 0.10 Mj^p and a radius of 1.79 ± 0.09 Rjup, 
WASP- 12b is the planet with the second largest radius reported to date. It is also the most 
heavily irradiated planet known, with an incident stellar flux at the substellar point of over 
9x10^ erg cm~^ s~^. In addition, model fits to its observed radial velocity and transit light 
curves suggest that the orbit of WASP-12b is slightly eccentric. All these attributes make 
WASP-12b one of the best targets to test current irradiated atmosphere and tidal heating 
models for extrasolar planets. 

In irradiated atmosphere model studies WASP- 12b is an extreme case even in the cat- 
egory of highly irradiated gas giants. Such highly irra diated planets are ex pected to show 
thermal inversions in their upper atmospheric layers (IBurrows et al.l 120081 ). although the 
chemicals responsible for such inversions remain unkno wn. TiO and VO molecules, which 



can a ct as strong optica l abso rbers, have been proposed ( IHubeny et al.ll2003l : iFortney et al. 



20081 ). but iDesert et al.l (120081 ) claim that the concentr ation of those molec ules in planetary 
atmospheres is too low (< 10~^ — 10"^ times solar). ISpiegel et al.l (120091 ) argue that TiO 
needs to be at least half the solar abundance to cause thermal inversions, and very high 
levels of macroscopic mixing are required to keep enough TiO in the upper atmosphere of 
planets. S2, and HS compounds have als o recently been sugg ested and then questioned 



as causes of the observed thermal inversions JZahnle et al.ll2009al lbh. 



I n the case of tidal heating, detailed ra od els are now being develo ped (e.g. iBodenheimer et al. 



200.1 iMiller et al.l 120091 : llbgui et al.l l2009l Jbl: llbgui fc Burrowsl l2009l ) to explain the inflated 
radius phenomenon observed in hot Jupiters, of which WASP-12b, with a radius over 40% 
larger than predicted by standard models, is also an extreme case. All models assume that 
the planetary orbits are slightly eccentric, and directly measuring those eccentricities is key 
not only to test the model hypotheses, but also to o btain infornaation about the planets' 
core mass and energy dissipation mechanisms (see e.g. llbgui et al.l 120091 ). 



Here we present the detection of the eclipse of WASP-12b in the 2;'-band (0.9 /xm), which 
gives the first measurement of the atmospheric emission of this planet, and the first direct 
estimation of its orbital eccentricity. Section 2 summarizes the observations and analysis of 
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the data. In Section 3 we compare the emission of the planet to models. The results are 
discussed in Section 4. 



2. Observations and Analyses 

We monitored WASP-12 [RA(J2000) =06:30:32.794, Dec(J2000)=+29:40:20.29, V=11.7] 
during two echpses, and under photometric conditions, on February 19 and October 18 2009 
UT. An additional attempt on October 30 2009 UT was lost due to weather. The data were 
collected with the SPICam instrument on the ARC's 3.5-meter telescope at Apache Point 
Observatory, using a SDSS z' filter with an effective central wavelength of ~0.9/im. SPICam 
is a backside-illuminated SlTc TK2048E 2048x2048 pixel CCD with 24 micron pixels, giving 
an unbinncd plate scale of 0.14 arc seconds per pixel and a field of view of 4.78 arc minutes 
square. The detector, cosmetically excellent and linear through the full A/D converter range, 
was binned 2x2, which gives a gain of 3.35 e~/ADU, a read noise of 1.9 DN/pixel, and a 48 
second read time. 

On February 19 wc monitored WASP-12 from 3:00 to 3:28 UT and from 3:54 to 7:10 
UT, losing coverage between 3:28 and 3:54 UT when the star reached a local altitude greater 
than 85°, the soft hmit of the telescope at that time. These observations yielded 1.20 hours 
of out-of-echpse and 2.45 hours of in-echpse coverage, at airmasses between 1.005-1.412. On 
October 18 we extended the altitude soft hmit of the telescope to 87° and covered the entire 
eclipse from 7:05 to 12:45 UT, yielding 2.73 hours of out-of-eclipse and 2.93 hours of in-eclipse 
coverage, with airmasses between 1.001-1.801. In both nights we defocused the telescope to 
a FWHM of ~ 2" to reduce pixel sensitivity variation effects, and also to allow for longer 
integration times, which minimized scintillation noise and optimized the duty cycle of the 
observations. Pointing changed by less than (x,y) = (4,7) pixels in the October 18 dataset, 
and by less than (x,y) = (3,12) pixels on February 19, with the images for this second night 
suffering a small gradual drift in the y direction throughout the night. Integration times 
ranged from 10 to 20 seconds. Taking into account Poisson, readout, and scintillation noise, 
the photometric precision on WASP-12 and other bright stars in the images ranged between 
0.07-0.15% per exposure on February 19, and between 0.05-0.09% per exposure on October 
18. 

The field of view of SPICam was centered at RA(J2000)=06:30:25, Dec(J2000)=+29:42:05 
and included WASP-12 and two other isolated stars at RA(J2000)=06:30:31.8, Dec(J2000)=+29:42:27 
and RA(J2000)=06:30:22.6, Dcc(J2000)=+29:44:42, with apparent brightness and B-V and 
J — K colors similar to the target. Each night's dataset was analyzed independently and the 
results combined in the end. The timing information was extracted from the headers of the 
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images and converted into Heliocentric Julian Days using the IRAF task setjd, which has 
been tested to provide sub-second timing accuracy. 

We corrected each image for bias-level and flatfield effects using standard IRAF routines. 
Dark current was neglegible. DAOPHOT-type aperture photometry was performed in each 
frame. We recorded the flux from the target and the comparison stars over a wide range of 
apertures and sky background annuli around each star. We used apertures between 2 and 
35 pixels in one-pixel steps during a first preliminary photometry pass, and 0.05 pixel steps 
in the final photometric extraction. To compute the sky background around each star we 
used variable width annuli, with inner radii between 35 and 60 pixels sampled in one-pixel 
steps. 

The best aperture and sky annuli combinations were selected by identifying the most 
stable (i.e. minimum standard deviation), differential light curves between each comparison 
and the target at phases out-of-eclipsqj. In the February 19 data, the best photometry 
results from an aperture radius of 14.7 pixels for both the target and the comparison stars, 
and sky annuli with a 52-pixel inner radius and 22-pixel wide. For the October 18 data, 17.9 
pixel apertures and sky annuli with a 45-pixel inner radius and 22-pixel wide produce the 
best photometry. 

The resultant differential light curves between the target and each comparison contain 
systematic trends that can be attributed to either atmospheric effects, such as airmass, see- 
ing, or sky brightness variations, or to instrumental effects, such as small changes in the 
location of the stars on the detector. Systematics can also be introduced by instrumental 
temperature or pressure changes, but those parameters are not monitored in SPICam. We 
modeled systematics for each light curve by fitting linear correlations between each param- 
eter (airmass, seeing, sky brightness variations, and target position) and the out-of-eclipse 
portions of the light curves. All detected trends are linear and there are no apparent residual 
color difference effects. The full light curves are then de-trended using those correlation fits. 
In the October 18 dataset, airmass effects are the dominant systematic, introducing a linear 
baseline trend with an amplitude in flux of 0.07%. The February 19 dataset also shows 
systematics with seeing and time with a total amplitude of also 0.07%. The systematics 
on this night were modeled using only the after-eclipse portion of the light curve, and we 
consider this dataset less reliable that the October 18 one. The 18 pre-ingress images col- 
lected between 3:00 and 3:38 UT suffer from a ~50 pixel position shift with respect to the 
rest of the images collected that night, which cannot be modeled using overall out-of-eclipse 



^We had to iterate on the out-of-echpse phase limits after finding that the echpse was centered at </) = 
0.51. Out-of-ecHpse was finally defined as phases 0<O.45 and (j)>0.57. 
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Fig. 1. — Top: De-trended differential light curves. Open and filled dots show, respectively, 
the Feb 19 and Oct 18 UT 2009 data. Bottom: Combined light curves binned by a factor 
of 12. Blue squares correspond to WASP- 12 and black dots to the differential light curve 
of the two comparison stars. The best fit models are shown as solid lines (for e = 0) and 
dashed lines (for e = 0.057). Both models produce the same depth and center phase, but 
the e = 0.057 model lasts 11.52 minutes longer. We attribute the flux jumps between phases 
0.475 and 0.5 to unremoved systematics. Notice that, although the systematics appear in 
both curves, the trends in each curve are not correlated in phase. 



systematics. We chose not to use those points in the flnal analysis. Correlations with the 
other parameters listed above are not signiflcant in any of the two datasets. 

Finally, we produce one light curve per night by combining the de-trended light curves 
of each comparison. The light curves are combined applying a weighted average based on 
the Poisson noise of the individual light curve points. The result is illustrated in Figure 1. 
The out-of-eclipse scatter of the combined light curves is 0.11% for the February 19 data 
and 0.09% for the October 18 data. De-trending signiflcantly improves the systematics, but 
some unidentifled residual noise sources remain, which we have not been able to fully model. 
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2.1. Eclipse detection and error estimation 



The two -night comb i ned lig ht curve contains 421 points between phases 0.413 and 0.596, 
based on the iHebb et al.l (120091 ) ephemerides. To establish the presence of the eclipse and its 
parameters, we fit the data to a grid of models generate d using the JASMINE code, which 
combines the iKippingl ( 120081 ) and iMandel fc Agoll ( 120021 ) algorithms to produce model light 
curves in the general case of eccentric orbits. The models do not include limb darkening 
and use as input parameters the orbital period, stellar and planetary radii, argument of the 
periastron, orbital inclination, stellar radial velocity amplitude and semi-major axis values 
derived by Hebb et al. (2009). The eccentricity is initially assumed to be e=0, which 
produces models with a total eclipse duration of 2.808 hours. The best fit model is found 
by minimization, with the depth, the central phase of the eclipse, and the out-of-eclipse 
differential flux as free parameters. 

First we fit the individual night light curves to ensure the eclipse signal is present in 
each dataset. The February 19 data give an eclipse depth of 0.100 ± 0.023%, while the 
derived echpse depth for the October 18 data is 0.068 ± 0.021%. The central phases are 
0=0.510 for the first eclipse and = 0.508 for the second. We assume the difference in depth 
is due to systematics we have not been able to properly model. The incomplete eclipse from 
February 19 might seem more prone to systematics, but our inspection of both datasets does 
not reveal stronger trends in that dataset. We therefore combined the data from both nights, 
weighting each light curve based on its out-of-eclipse scatter. 

The result of the combined light curve analysis is the detection of an eclipse with a 
depth of 0.082 ± 0.015% and centered at orbital phase = 0.51, as shown in Figure 1. The 
reduced of the fit is 0.952. The error in the eclipse depth is computed using the equation 
'^depth — "^w/^ + "^r) wherc aw is the scatter per out-of-eclipse dat a point and de scribes 



the red noise. The cr^ is estimated with the binning technique by iPont et al.l (120061 ) to be 
1.5x10"^ when binning on timescales up to the ingress and egress duration of about 20 
minutes. 

We investigate to what extent the uncertainties in the system's parameters affect our 
eclipse depth and central phase results. Varying the impa ct paramete r , plan et-to-star ratio, 
and scale of the system by la of the reported values in iHebb et al.l (120091 ) . the measured 
eclipse depth changes only by 0.004% or 0.27adepth, while the central phase remains un- 
changed. Our result is therefore largely independent of the adopted system parameters. 

We performed three other tests to confirm the echpse detection in a manner simi 



lar to previously r eported echpse results (IDeming et al.l l2005l : ISing fc Lopez-Moraled 12009 



Rogers et al.ll2009l ). From the average of the 125 out-of-eclipse light curve data points versus 
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the 228 in-eclipse points (only points where the planet is fully eclipsed, adopting {/)=0.51 as 
the central eclipse phase), we measure an eclipse depth of 0.080 ± 0.015%. We further check 
the detection by producing histograms of the normalized light curve flux distribution in the 
in-eclipse and out-of-eclipse portions of the light curve. The result, illustrated in Figure 2, 
shows how the flux distribution of in-eclipse points is shifted by 0.00082 with respect to the 
out-of-eclipse flux distribution, centered at zero. The results of the last test, where we use a 
new set eclipse of models with a duration corresponding to e = 0.057, fix the out-of-echpse 
baseline to zero, and leave both depth and central phase of the eclipse as a free parameters, 
are shown in Figure 3. The two top panels in the figure show 1-D cuts of each parameter 
through the z (normalized x^)) space, where 2=0 gives the best fit model and z=\ defines 
the 1 ( 7 confidence interval of the result (see definition of z in §3 of iLopez-Morales &: Bonanos 



20081 : iBehrl |2003[ ) . However, inspection of the eclipse depth in Figure 1 reveals that the la 



errors from this method (±0.036%) are too large. To estimate more realistic errors for the 
central phase, we generate contours plots (bottom panel in Fig. 3), where the la eclipse 
depth error has been fixed to the 0.015% value derived above. The resultant center phase 
= 0.5100^Q;oQg^. We also applied praye r-bead, bootstrapping, and Markov-C hain Monte 



IS 



Carlo (MCMC) error analysis techniques ( iGillon et al.l 120071 : ISouthworthI 120081 ) to estimate 
the errors in the central phase of the eclipse. The MCMC analysis contained 2.85x10^ links, 
and we adopted t he 1.5x10"^ photometric error (~1.6 times the formal error), given by the 



Pont et al.l (|2006[ ) binning technique. All these give central phases within = 0.5100 



+0.0030 
-0.0036- 

The larger error in the central phase given by the normalized method suggests the presence 
of correlated noise in the data, which the other methods might not be correctly accounting 
for. We therefore adopt these larger errorbars in our final estimation of the central phase. 



2.2. Eccentricity 



The eccentricity e o f WASP- 12b was c alculated from the measured central phase shift 
value using eq. (6) from IWallenquistI (Il950l ). 

TT (t2 - tl - P/2) 



ecoscj 



1) 



P 1 + csc^ i ' 

where P, i and u are, respectively, the orbital period, inclination, and periastron angle 
of the system, and ^2 ~ is the time difference betwee n transit and eclip se. In our case 
t2 — t\ = 0.51P. Using the values of P, i and u from iHebb et al.l (l2009l ). we derive an 
e = 0.057lo;o34, which agrees with the non-zero eccentricity result reported by these authors. 
This eccentricity can be in principle explained if 1) the system is too young to have already 
circularized, 2) there are additional bodies in the system p umping the eccentricity of WASP- 
12b, or 3) the tidal dissipation factor Qp (lGoldreichlll963l ) of WASP-12b is several orders of 
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magnitude larger than Jupiter's, estimated to be between 6 x 10^ and 2 x 10® flYoder &: Peak 



198l|). 



3. Comparison with atmospheric models 

We compare the observed 2;'-band flux of WASP- 12b to simple blackbody models and to 
expectedly more realistic radiative-convective models of irradiate d planetary at n aosph eres in 
chemical equilibrium, following the same procedure described in [Rogers et al.l (120091 ) . The 
results are shown in Figures 4 and 5. 



In the simplistic blackbody approximation, a 0.082±0.015% deep eclipse corresponds 
to a z'-band brightness temperature of T^/ ~ 3028^^^o-^5 slightly lower than the planet's 
equilibrium temperature o f Tp ^ 3129.ft" assuming zero Bo nd albedo {Ab = 0) and no energy 



reradiation (/ = 2/3) (see iLopez- Morales fc Seagerl 120071 ) . However, when the thermal and 



reflected flux of the planet are included, different combinations oi Ab and / can yield the 
same eclipse depth, as illustrated in Figure 4. From that flgure we can constrain the energy 
redistribution factor to / > 0.585 ±0.080, but the albedo is not well constrained. Assuming 
a maximum Ab < 0.4, the temperature of the day-side of WASP-12b is Tp > 

The more realistic atmospheric models are derived from self-consistent coupled ra- 
diative transfer and chemic al equilibrium c a lculat ions, based on the models described in 
Sudarskv et al.l JgOOol. l2003h. iHubenv et all J2003h and burrows et aP J2005I . l20od . l2008h 



Rogers et al.ll2009l . for details). We generate models with and without thermal inversion 



{see 

layers, by adding an unidenfltied optical absorber between 0.43 and 1.0 /im, with different 
level of opacity k'. The opacity of the absorber varies parabolically with frequency, with a 
peak value of = 0.25 cm^ g^~^- As Figure 5 shows, models with and without extra ab- 
sorbers produce similar flts to the observed z'-band flux. The best model without absorber 
has a Pn = 0.30. The best model with an extra absorber has a -P„ = 0.1 and = 0.1 cm^ 
gr~^. Observations at other wavelengths are necessary to further constrain the models. 



^Pn ~ and Pn = 0.5 correspond, respectively, to / = 2/3 and f = 1/4, however there is not well defined 
Pn — f relation for intermediate values since the physical models account for atmospheric parameters (e.g. 



pressure, opacity) in a way different than blackbody models. 
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Discussion and Conclusions 



This first detection of tiie eclip s e of W ASP-12b agrees witli tlie sliglit eccentricity of tlie 
planet's orbit found by lHebb et al.l (120091 ) . and places initial constraints to its atmospheric 
characteristics. 

The presence of other bodies in the system can be tested via radial veloc i ty or transit 
timing variation observations, although the current RV curve by iHebb et al.l (120091 ) shows 
no evidence of additional planets, unless they are in very long orbits. 

One would expect that if extra absorbers are present in the upper atmosphere of the 
planet in gaseous form, they might give rise to thermal inversion layers. However, as Figure 
5 illustrates, the observed 0.9 /xm eclipse depth can be fit equally well by a model without 
extra absorbers. Additional observations at longer wavelengths, specially longer than ~ 4.0 
fim, will break that model degeneracy. Observations at wavelengths below ~ 0.6 fim will 
also place better constraints on As- 

Note: While this paper was in the final revision stages , two new secondary eclipse 
observations from Spitzer were reported by ICampo et al.l (120101 ) . They find secondary eclipse 
central phases consistent with a circular orbit. We have carefully reviewed our analysis and 
still arrive to a slightly eccentric orbit, although the result is only s i gnifica nt at the 1.4-1.6(7 
confidence level. Possible explanations proposed by ICampo et al.l (120101 ) for the measured 
eccentricity difference are orbital precession or wavelength-dependent brightness variations 
across the surface of the planet that would s hift t he center of the eclipse. A similar effect 
has been recently suggested by ISwain et al.l (120101 ) on the surface of HD189733b. Further 
observations are needed to establish whether this discrepancy is a data artifact or a real 
effect. 
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Fig. 2. — Normalized flux histograms of the in-eclipse (dotted red line) and out-of-eclipse 
(solid line) portions of the WASP-12 light curve in Figure 1. The bin width is 0.00082 in 
differential flux, coincident with the detected eclipse depth. 



0.00 0.04 0.07 0.11 0.14 0.46 0.48 0.50 0.52 0.54 0.56 

Eclipse Depth (%) Center Phase 




0.48 0.50 0.52 0.54 

Center phase 



Fig. 3. — Top: Model eclipse 1-D cuts through the normalized-x^ parameter space for the 
eclipse depth and center phase. z=Q gives the best fit and z=l shows the la confidence 
interval. Bottom: Contour plot of the eclipse depth versus the center phase, where the la 
depth error has been fixed to 0.015%. The best model fit is indicated by the cross symbol at 
0=0.51 and depth = 0.082%, together with la to 3cr confidence contours. The star symbol 
at 0=0.50 corresponds to a circular orbit. 
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Fig. 4. — Values oi Ab and / that reproduce the observed 2;'-band eclipse depth of WASP- 
12b, assuming the planet emits as a blackbody. The shaded area highlights the region of 
allowed / values (1/4 — 2/3). The short and long dashed lines delimit, respectively, the la 
and 2a confidence regions of the result. 
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Fig. 5. — Comparison of the eclipse depth (planet-to-star flux ratio) to models. The green 
line shows the As = 0.4, / = 2/3 blackbody model from Figure 4. The blue and red lines 
show, respectively, the best fit model for an atmosphere with no extra absorber, and with 
an extra absorber of opacity k' between 0.43 and 1.0 /im. The black thin line at the bottom 
indicates the SPICam plus SDSS 2;'-band filter response. See §3 for more details. 



